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We show that hghtly doped YBa2Cu306+i/ has muhiple holon pockets with different areas, which 
lead to multiple frequencies of magnetic quantum oscillations. Using neutron scattering data on 
incommensurate spin ordering we determine these areas, which yields frequencies in good agreement 
with experiments. Divergence of the effective mass observed in magnetic quantum oscillations 
indicates a quantum phase transition at the oxygen content y « 0.48. We argue that the transition 
is the onset of quasistatic incommensurate magnetic order predicted by theory and observed in 
neutron scattering. 

PACS numbers: 74.25.Jb 71.18.+y 74.72. Gh 75.30.Fv 



Topology of the Fermi surfaee is one of the central 
problems in the physics of cuprate superconductors. An 
undopcd cuprate is antiferromagnetic Mott insulator. At 
a sufficiently small doping the system can be described by 
the extended two-dimensional t-J model [J, Q that pre- 
dicts small holon pockets located at the "nodal points" . 
We use the term "holon" to stress that the quasiparti- 
cle does not carry usual spin, with more details to be 
explained later. On the other hand there is no doubt 
that at a sufficiently large doping the system behaves like 
a normal metal with a large Fermi surface and normal 
quasiparticles which carry simultaneously electric charge 
and spin 5=1/2. Description of the transition between 
these two regimes is an open theoretical problem. 

Angle resolved photoemission spectroscopy (ARPES) 
studies indicate a large Fermi surface in overdoped regime 
and formation of Fermi arcs in the underdoped regime [H . 
On the other hand, recent magnetic quantum oscillation 
(MQO) data taken in YBa2Cu306+y(YBCO) in- 

dicate small Fermi pockets at y < 0.66 that corresponds 
to doping X < 0.125. Throughout the paper we denote 
doping by x and use results of Ref.ll3|to relate the oxygen 
content y with the doping level x. Very recent ARPES 
studies IJ, |l5| also give some indications of small Fermi 
pockets. A theoretical interpretation of ARPES in terms 
of holons is a fairly involved issue. In a photoemission 
measurement partially separated spin and charge have 
to recombine to form a physical electron. So far the re- 
combination amplitude has only been calculated for the 
parent Mott insulator In contrast, MQO is sensi- 

tive mainly to the electric charge, hence interpretation of 
MQO in terms of holons is straightforward. With tilted 
magnetic field MQO can also probe spin of the quasipar- 
ticle. A very recent measurement of this kind [lO| indi- 
cates suppression of spin and hence supports the holon 
picture. 

The sign of the Hall coefficient measured in Ref. i5i cor- 
responds to electron pockets instead of hole ones. How- 
ever, the field applied in the experiments is smaller than 



the critical field Hc2, hence the contribution from the 
vortex liquid has to be taken into account, which may 
alter sign of the Hall coefficient [l7, 18|. In the present 
work we consider only hole pockets. 

In this paper we address the following issues related to 
MQO: (l)Typically more than one MQO frequency is ob- 
served. What is the origin for the multiple frequencies? 

(2) The main MQO frequency corresponds to the area of 
the pocket about 1.8% of the Brillouin zone. How to rec- 
oncile this very small area with the Luttinger sum rule? 

(3) Divergence of the effective mass has been observed [11 1 
at doping about 9%. What is the physical origin for this 
quantum critical point (QCP)? 

We will relate MQO with recent neutron scattering ob- 
servations }19l - i23j . These data demonstrate incommen- 
surate spin ordering pined to the tetragonal a* direc- 
tion. Depending on doping the ordering can be static or 
dynamic. The QCP separating regions of dynamic and 
static ordering is located at xqcp ~ 0.09. The spin-wave 
pseudogap Agw is opened at x > xqcp- The data [l^ 
on YBa2Cu306.5 {x « 0.1) show Agw ~ 10 meV, while 
in YBaaCuaOe.e {x ~ 0.12) A^w 20 meV [Ifl. On 
the other hand the quasistatic scattering has been ob- 
served at x < xqcp where Agw = 0. In YBa2Cu306.45 
(x « 0.085) the signal is very weak [2l|, and it is 
much bigger ^ in YBaaCugOg.as {x w 0.065). 

MQO in a magnetic field B is described as cos(27r^ + 
(j>). The period F and the area enclosed b y a trajectory 
in the momentum space Ak are related as [2J] 
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The typical measured value of F in YBCO is about 500 T, 
which gives a ratio of A^ to the total area of the Brillouin 
zone A^i = (27r/a)^ (a = 3.81 A is the lattice spacing) 
to be Ak/Anz ~ 0.0176. 

Our analysis of MQO is based on the spin spiral the- 
ory of a lightly doped Mott insulator. While the idea 
of the spin spiral has been suggested quite some time 



2 



ago [25 -28 1, the consistent theory with account of quan- 
tum fluctuations has been developed only recently pof . 
Generalization of the spiral theory to bilayer YBCO has 
been discussed in Ref. [s^. The theory is based on ex- 
pansion in powers of doping a;, so it is parametrically 
justified at a: <C 1. The theory is formulated in terms 
of the bosonic n- field (n^ = 1) that describes the stag- 
gered component of the copper spins, and in terms of 
fermionic holons ip. The holon field carries electric 
charge and it has a pseudospin that originates from two 
sublatticcs. Minima of the holon dispersion are at the 
nodal points qo = (±7r/2, ±7r/2), so there arc holons 
of two types ipa^ a = 1,2, corresponding to two nodal 
directions. The dispersion in a pocket is somewhat 
anisotropic, but for simplicity wc use here the isotropic 
approximation, e (p) « /3p^/2, where p = q— qo. Lattice 
spacing of the square lattice is set to be equal to unity 
a = 3.8IA— >1. The effective Lagrangian for a single 
layer system reads [2^ 

a 

- ^ie{V)^^ + V2g{iPia^^) ■ [n x (e„ • V)r?]} , 
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First two terms in the Lagrangian represent the usual 
nonlinear a model, with magnetic susceptibility XJ- ~ 
0.066J and spin stiffness ps ~ 0.175J, where J « 
130 meV is the antiferromagnetic exchange in the par- 
ent Mott insulator. The extended t-J model predicts the 
following values of the coupling constant and the inverse 
effective mass, g ss l.OJ, /3 sa 2.5J (m* = LBrrie). Below 
wc will not use this value of /?, instead we will treat 13 
as a fitting parameter. The pseudospin operator is ^tr, 
and = (1/^2, ±1/^2) is a unit vector orthogonal to 
the face of the magnetic Brillouin zone (MBZ) where the 
holon is located. Note that usage of the MBZ notations 
does not imply that there is a long range magnetic or- 
der. This is just a convenient way to avoid double count- 
ing of degrees of freedom, and the pseudospin accounts 
for doubling of the area of MBZ. Due to the MBZ no- 
tations one should consider two full pockets located at 
qo = (-l-7r/2, — 7r/2) and qo = (+7r/2, +7r/2), hence the 
index a takes two values. We stress that it does not 
matter if the ground state expectation value of the n- 
field is nonzero, (n) ^ (magnetic ordering), or zero, 
(n) =0. The only condition for validity of Eq.([2]) is that 
dynamic fluctuations of the n-field are sufficiently slow. 
Typical energy scale of the n-field quantum fluctuations 
is i?cross fx x^^"^ (position of the neck of the "hour glass" 
spin wave dispersion), as discussed in Ref. [2^ and it must 
be small compared to the holon Fermi energy oc x. 
The inequality -Ecross <C Ci? is valid up to optimal doping. 



X « 0.15, below which Eq.([2]) is parametrically justified. 

To account for the interaction with magnetic field, we 
make the following modifications in Eq.^ [2^: A mag- 
netic vector potential A is included in the long deriva- 
tive V ^ V — |A, kinetic energy of n— field is modi- 
fied n — >■ — [ri X B]^ , and an extra term SCb = 

^ (B ■ n)'4'l, ■ ^) V'a is included. Here we include the Bohr 
magneton in definition of the magnetic field, 2p,-B,B B. 
The precise meaning of spin-charge separation is clear 
under this context: In a normal Fermi liquid the spin in- 
teraction is 5C^^ = ^ipl^B ■ aipa- In contrast, for a spin- 
spiral state the g-term in Eq. ^ enforces ri _L ct, hence the 
expectation value of 6£b goes to zero, {{B ■n){a -n)) = 0. 
Thus the pseudospin does not interact with magnetic 
field in the first order in B and this is the meaning of 
the spinless nature of the holon. 

In the effective action ^ the p^-term. the e(7')-term 
and the g-tcrm are important in semiclassical approxi- 
mation. All other terms contain time derivatives, which 
are important only for quantum fiuctuations and we dis- 
regard them for now. For the single layer case, we choose 
the plane of the coplanar spin spiral to be the a;j/-plane 
in the spin space, n ^ (cosQ • r,sinQ • r,0). Due to the 
g-term in Eq.(l2|) the psedospin of a holon is quantized 
along z-axis. Energy of the holon is e(p) = zLgQ + /3p'^/2, 
where ±gQ is the pseudospin splitting. The splitting be- 
tween two pseudospin branches is shown in the left panel 
of FigHJ Minimization of the semiclassical energy shows 




FIG. 1: (color online) Filling of holon bands. Left panel: sin- 
gle layer La2-iSr2,Cu04 (LSCO). Three right panels: double 
layer YBCO in three different regimes. The solid and the 
dashed line correspond to different pseudospin projections, 
the splitting is ±gQ. In the YBCO panels left/right parts 
show bonding/antibonding bands, the splitting is ±A/2. 

that Q is directed along the CuO bond [Q cx (1,0) or 
Q (X (0, 1)] and has value Q = j-x. The upper pseu- 
dospin branch of the holon dispersion is always empty, as 
shown in the left panel, Fig[Tl hence the area of the holon 
pocket is = ^x. Quantum fiuctuations reduce the 
static value of spin, |(n)| < 1, moreover, at a; > xqcp ~ 
0.1 the expectation value of ft vanishes, (ft) = 0, and the 
spiral becomes fully dynamic (29j . However, this does 
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not influence the value of the incommensurate vector Q 
because the fluctuations are slow. The point is that the 
semiclassical analysis is based on relatively short distance 
and time, with corresponding typical momenta and en- 
ergies Q < (Jsomi < PF OC v^, Wsemi ~ EF OC X. Quautum 

fluctuations come from smaller momentum/energy scales 
Q'quant < Q, ^^quant oc x^/^, as discussed in Ref. |2^. 

In the double layer case, one has to include the in- 
ter layer hopping in the effective action ([2]). As a re- 
sult, holon wave function ij; attains bonding/antibonding 
index with respect to the interlayer hybridization (soj . 
The holon energy at each of the two nodal points reads 
e(p) = /?p^/2 ± gQ ± A/2, where ±gQ is the pseudospin 
splitting and ±A/2 is the antibonding(a)/bonding(b) 
splitting. Effectively there are four different bands per 
nodal direction (&— , b+, a—, a-|-), altogether eight bands. 
The filling configuration of these bands is determined by 
minimizing the semiclassical energy, which yields three 
different doping regimes [sof 



Q = 



0.15 



1) .T < Xo , 

2) xo < X < xi , Q ~ 

3) X > xi , 



g X- A/inp) 
Ps 3 - 2A 

Q ^ gx/ps , 



(3) 



where A ~ 2g'^ / {t: (3 p s) . The points xq and xi are Lifshitz 
points. Filling configuration of these three regimes are 
shown in Fig[TJ Areas of the filled holon pockets are 



1) A6-/Abz 

2) A6_Mbz 

3) A_Mbz 

^a-MsZ 



Ab+/ABz = .t/2 
.x/3 + .gO/(37r/?)4 
x/3-2.gQ/(37r/3) 
x/3 + .gQ/(37r/3)- 
x/2-f A/(47r/3) 
x/2- A/(47r/3) . 



A/(67r/3) 
f A/67r/3) 
A/(37r/3) 



(4) 



Note that the total occupied area is 2x because there are 
two layers. From ^ we conclude that there is one MQO 
frequency in the first regime, three frequencies in the sec- 
ond regime, and two frequencies in the third regime. Let 
us consider two scenarios for the doping dependence of 
the hybridization splitting: (A) a constant gap A = Aq, 
(B) A = Ao(l -I- vx), since the splitting A is due to 



a linear de- 



tunneling via interlayer oxygen chains [31 1 
pendence on doping is possible. 

Values of the incommensurate wave vector Q deter- 
mined by neutron scattering 3"21| are presented in 
FigEl Comparing FigOwith Eqs.dS]) one finds xq ~ 0.06, 
xi « 0.13. Using Eqs.(I5]), we can determine parameters 
of the model. We fix g = J and ps ~ 0.175J, as they 
are predicted by the extended t-J model. Within the 
scenario (A) the fit gives, A = 1.23, /3 = 2.95J (m* = 
1.35me), Aq = 0.556J. Within the scenario (B), the data 
from Fig[2] arc not sufficient to determine the additional 
parameter v. However, we can assume that value of A in 
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FIG. 2: (color online) Incommensurate wave vector versus 
doping. The blue square [3, the red circle [l^, and the red 
diamond [U] show neutron scattering data. The solid line 
shows fit of the data using Eqs.((3]|. 

YBCO and LSCO is the same, and use the value A — 1.31 
obtained in Ref. [2^ from fitting the data for LSCO. This 
gives 13 = 2.78J (m* = 1.43me), Aq = 0.37J, v = 6.9 
within the scenario (B). Note that values of (3 in both fits 
are close to that predicted by the extended t-J model, as 
demonstrated in the paragraph after Eq. ([2]) , and the cor- 
responding effective masses are close to that measured in 
MQO away from the QCP [ul [I!], m* = 1.6 ± O.lrrie. 
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FIG. 3: (color online) Frequencies of MQO in Tesla versus 
doping. Solid lines show the theoretical prediction within 
the scenario (A), dashed lines show the prediction within the 
scenario (B). Predictions (A) and (B) coincide in the lower 
part of the figure. Experimental data are shown by symbols: 
Ref. 0] - blue circles, Ref. - black triangles down, Ref. 
- black triangles up. Full symbols correspond to maximum 
intensity lines. 

Frequencies of MQO are determined by these param- 
eters via Eqs.dl]) and ([T]). The results are presented in 
FigEl Difference between the scenario (A) (solid lines) 
and the scenario (B) (dashed lines) is fairly small. There 
is quite a reasonable agreement with experimental data, 
especially having in mind that the theory has no fitting 
parameters related to MQO. Even a high frequency point 
« 1690r is reproduced. Notice that MQO always con- 
tain higher harmonics. Most likely the blue circle at 
w HOOT is such a harmonic [ll|. Filling of each pocket 
is only a fraction of the total doping, this explains why 
the filling extracted from one single frequency seems to 
violate the Luttinger sum rule. The main frequency at 
X sa 0.09 comes from the b+ pocket, which has Fermi 



4 



energy w 40 meV, and the number of filled Landau levels 
at B = 60Tesla is N = ep/fv^c ~ 6, ujc being the syn- 
chrotron frequency. Amplitude of MQO is proportional 
to exp {— TrF/ZzWc}, Ref. (s^ . Assuming that the expo- 
nent is ~ 1, and relating the impurity broadening to the 
mean free path and to the Fermi velocity, F = hvp/l, we 
estimate the holon mean free path I ^ 70a ^ 270 A. 

A divergence of the effective mass at approaching 
y = 0.49 from higher doping has been reported The 
effective mass has been extracted via temperature de- 
pendence of the MQO signal. The divergence was inter- 
preted [Uj as a QCP due to a metal-insulator transition. 
Here we suggest an alternative explanation for this QCP, 
based on the current theory and on the evidences from 
neutron scattering. Besides the Lifshitz points xq and 
xi, the spiral theory also predicts a QCP at xqcp ~ 0.1 
which separates regions of static and dynamic spin spi- 



ral. The neutron scattering data 1^22| demonstrate 
this QCP in YBCO at doping level xqcp ~ 0.09, see 
discussion in the introduction. As one should expect, 
the precise position of this QCP depends on the applied 
magnetic field. The quasistatic neutron scattering in 
YBa2Cu3 06.45 is enhanced in the field p^ . which indi- 
cates that magnetic field shifts the QCP towards higher 
doping. Most likely in the field 60-80 Tesla used in 
MQO; the QCP is located between YBa2Cu306.47 and 
YBa2Cu3 06.49. The point is that the field is hardly suf- 
ficient to close the spin-wave pseudogap Agw ~ 10 meV 
in YBasCuaOg.s [li- The QCP is driven by the long- 
wave-length quantum fiuctuations, while holon pockets 
are formed at a shorter semiclassical scale. Therefore 
MQO frequencies are not sensitive to the magnetic QCP. 

Effective mass in a quantum field theory always de- 
pends on the momentum transfer q, m* — ^ m*. The 
effective mass extracted from fitting neutron scattering 
data is relevant to the semiclassical scale q ^ pp (x. ^/x. 
This mass does not "see" the magnetic QCP. However, 
the amplitude of the MQO signal is formed at the length 
scale I ~ 70a ~ 270 A, which corresponds to a very small 
q and the corresponding effective mass is sensitive to the 
magnetic QCP. It is known [s^l that the quasiparticle 
residue Z vanishes at a magnetic QCP. Since m* cx 1/Z, 
this implies that the effective mass is diverging. Alter- 
natively one can say that the holon scattering rate from 
critical magnons is diverging. This gives a natural expla- 
nation of the divergence observed in Ref. [ll|, and also 



explains why position of the QCP observed in neutron 
scattering coincides with that observed in MQO: they 
are the same QCP. 

In summary, we explained the MQO frequencies based 
on the spin-spiral theory of lightly doped Mott insulator. 
Fit of the incommensurate neutron scattering data shown 
in Fig 12] determines free parameters of the theory, and al- 
lows us to calculate frequencies of MQO as functions of 
doping. In the doping region 0.06 < a; < 0.13 the theory 
predicts three different frequencies, which yields a fairly 



good agreement with current experimental data, Fig|3l 
and the Luttinger sum rule is reconciled. We also ar- 
gue that the quantum critical point observed in magnetic 
quantum oscillations (divergence of the effective mass) is 
the same quantum critical point that is observed in neu- 
tron scattering, which signatures the onset of the static 
incommensurate magnetic order. 

We thank D. Haug, V. Hinkov and B. Keimer for im- 
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data. We acknowledge useful comments by S. Borisenko 
and C. Zhang. 
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